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We describe a model for the thermodynamics and dynamics of glass-forming liquids in terms 
of excitations from an ideal glass state to a Gaussian manifold of configurationally excited states. 
The quantitative fit of this three parameter model to the experimental data on excess entropy 
and heat capacity shows that "fragile" behavior, indicated by a sharply rising excess heat capacity 
as the glass transition is approached from above, occurs in anticipation of a first-order transition 
- usually hidden below the glass transition - to a "strong" liquid state of low excess entropy. 
The distinction between fragile and strong behavior of glass-formers is traced back to an order of 
magnitude difference in the Gaussian width of their excitation energies. Simple relations connect the 
excess heat capacity to the Gaussian width parameter, and the liquid-liquid transition temperature, 
and strong, testable, predictions concerning the distinct properties of energy landscape for fragile 
liquids are made. The dynamic model relates relaxation to a hierarchical sequence of excitation 
events each involving the probability of accumulating sufficient kinetic energy on a separate excitable 
unit. Super- Arrhenius behavior of the relaxation rates, and the known correlation of kinetic with 
thermodynamic fragility, both follow from the way the rugged landscape induces fluctuations in 
the partitioning of energy between vibrational and configurational manifolds. A relation is derived 
in which the configurational heat capacity, rather than the configurational entropy of the Adam 
Gibbs equation, controls the temperature dependence of the relaxation times, and this gives a 
comparable account of the experimental observations. The familiar coincidence of zero mobility and 
Kauzmann temperatures is obtained as an approximate extrapolation of the theoretical equations. 
The comparison of the fits to excess thermodynamic properties of laboratory glass-formers, and to 
configurational thermodynamics from simulations, reveals that the major portion of the excitation 
entropy responsible for fragile behavior resides in the low-frequency vibrational density of states. 
The thermodynamic transition predicted for fragile liquids emerges from beneath the glass transition 
in case of laboratory water, and the unusual heat capacity behavior observed for this much studied 
liquid can be closely reproduced by the model. 



I. INTRODUCTION 

Viscous liquids close to the glass transition are usu- 
ally characterized by a broad distributions of (their) 
relaxation times (stretched exponential kinetics) and 
stronger than Arrhenius (super- Arrhenius) dependence 
of the relaxation times on temperature pi£ The stretched 
exponential kinetics is often represented by Kohlrausch- 
Williams- Watts (KWW) relaxation function 



,KWW 



(i) = exp [-(t/rf] 



(1) 



where (3 < 1 is a stretching exponent and </> KWW (i) is a 
normalized function representing some relaxing property. 
The temperature dependence of the relaxation time r in 
this equation has been given many forms ^ but is most 
commonly represented by the empirical Vogel-Fulcher- 
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Tammann (VFT) law 

ln(r/ro)=i)To/(T-To), 



(2) 
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were tq is a characteristic time for liquid quasi-latticc 
vibrations. According to this equation, the relaxation 
time diverges at the VFT temperature To below the glass 
transition temperature T g ; Tq is close to T g for fragile 
liquids with super-Arrhenius kinetics and tends to zero 
for strong liquids with nearly Arrhenius relaxation. 

At sufficiently low temperature, below the range of 
temperatures described by the mode-coupling theory,— 
a liquid can be characterized by a set of minima (inher- 
ent structures) of the potential energy landscape divided 
by Stillinger into "basins of attraction" £ The relaxation 
then occurs by a sequence of activated transitions^ be- 
tween basins or collections of basins separated by low 
barriers (metabasins)^^^ While opinions of its origin 
differ, there is much evidence for a dynamic crossover 
temperature in fragile liquids, usually corresponding to 
the critical temperature of mode coupling theory, where 
relaxation times are about 10~ 7 sJ^ Phenomenological 
theories, usually operating below this critical tempera- 
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ture, differ in their attribution of the driving force re- 
sponsible for activated events. 

The Adam-Gibbs theory^ 1 - suggests that entropy con- 
trols the activation process. According to this line of 
thought ) 12 ' 13 ' 14 the divergence of the relaxation time at 
To is caused by a decrease of the number of available 
states related to the configurational entropy S C (T). The 
relaxation time is given by the Adam-Gibbs relation 



ln(r/r ) = 



A 



TS C (T) 



(3) 



which anticipates a growing length-scale £ attributed to 
cooperatively rearranging regions. £ diverges as 5 C (T) _1 
on approach to the Kauzmann temperature Tk at which 
the configurational entropy vanishes ! 11 ! 15 ' 16 



S C {T K ) = 0. 



(4) 



The free volume theor} r 17 i 18 also considers entropy as 
the driving force for the glass transition in terms of the 
volume available for reorganizing the liquid. This "free" 
volume becomes increasingly scarce on cooling, eventu- 
ally leading to a divergent relaxation time at the VFT 
temperature. The analysis of T, P viscosity data, how- 
ever, indicates that it is temperature and not density 
that is the primary factor behind the super- Arrhcnius 
kinetics J£ 

The coefficient A in Eq. © can be related to system 
properties within the concept of entropic droplet i 20 ' 21 ' 22 
Following arguments by Bouchaud and Biroli^ 3 - and by 
Lubchcnko and Wolynes,— the probability of creation of 
a droplet of size £ is a competition of the surface, oc £ 2 , 
and entropic bulk, oc £ 3 , effects: 



p(o ex ex P (~a(oe + s c e 



(5) 



The probability minimizes at the stationary point of the 
exponent. Scaling of the surface tension with the droplet 
size of the form er(£) oc £ -1 / 2 results in the Adam-Gibbs 
law for the relaxation time, which, in this concept, cor- 
responds to the time of creating a mobile region rich 
in configurational entropy. This picture, however, does 
not directly address the question of the origin of viscous 
fiowp assuming that the mosaic structure of the dynam- 
ically exchanging mobile and immobile regions will have 
all the properties necessary to facilitate shear relaxation. 
The divergent length-scale of mosaic regions scales as^ 5 - 
£ oc S C (T)~ 2 / 3 in contrast to 5 C (T) _1 scaling in the 
Adam-Gibbs theory. 

The Adam-Gibbs relation can be tested by 
calorimetr y 26 ' 27 or by using relaxation times 
and configurational entropies from computer 
experiments.— i22i22i2ii22i22iM In the former case, 
one has to assume that S c can be approximated by 
the excess entropy of a liquid over the crystal S cx , 
although this is known to be a poor approximation in 
many cases ^ 5 -^ 6 - ' 37 ' 38 ^ Nevertheless, the Adam-Gibbs 
equation is known to work well for both 5 ex and S c , 
although it does not produce perfectly straight lines for 



ln(r) vs l/TS cyi , in particular for fragile liquids. The 
use of S CK in the Adam-Gibbs equation generally fails 
above the temperature Tb at which the dynamics change 
character (e.g. bifurcation into a and (3 processes)^ 
The Adam-Gibbs equation with S' cx used for S c is 
often equivalent to the VFT equation because many 
glass- former o 27 ' 40 empirically follow the \/T entropy 
decay: 



S CK {T) = S (1 - T™ p /T) 



(6) 



where T^ p is what we will call the experimental Kauz- 
mann temperature obtained by extrapolating S CX ~(T) to 
zero, in contrast to the thermodynamic Kauzmann tem- 
perature Tk defined by Eq. Q. Equation © often ap- 
plies to configurational entropies obtained from simula- 
tions, although bilinear in l/T forms have also been used 
to fit the data4^ 

A strong experimental argument in favor of the 
Adams-Gibbs picture is the near equality^ 3 - of the Kauz- 
mann temperature T K and the VFT temperature Tq. 
Further evidence is provided by the good account it 
gives of the pressure dependence of the glass transition 
temperature. 44 ' 45 ' 46 However, the theoretical relevance of 
the concept of nucleation of cooperatively rearranging re- 
gions is not clear. Given the fact that laboratory glass- 
formers and model fluids equilibrated at higher temper- 
atures in computer experiment approximately follow Eq. 
([B|, one wonders to what extent the combination of Eqs. 
([3]) and ([6]) is just a successful mathematical relation, 
reproducing the VFT law in Eq. @, and whether acti- 
vated events in supercooled liquids are really driven by 
the entropy. The resolution of this problem is important 
from a general perspective since the Adam-Gibbs picture 
puts supercooled liquids in a unique position within the 
more general problem of activated events in condensed 
matter. For most cases, including the vast majority of 
chemical reactions, excess kinetic energy at a transform- 
ing unit or a molecular mode, and not the entropy, is 
the driving force which lifts the system to the top of the 
activation barrier. This, more traditional, view of relax- 
ation of supercooled liquids is advocated by models that 
consider kinetic energy or enthalpy as the driving force 
of activated transitions^ 

In models of activation controlled by the kinetic en- 
ergy, one considers excitations to a common high-energy 
level above the "top" of the landscape corresponding to 
high-temperature diffusion^ 8 - The energy of the start- 
ing point at a basin minimum is treated as a ran- 
dom variable within trap models^^ or random-walk 
models . 50 ' 51 ' 52 All such models result in activated ki- 
netics with a temperature-dependent activation energy. 
However, since the probability of transition is finite at 
each temperature, no divergent relaxation time appears 
in these models. 

A potential advantage of the energy models is the 
opportunity to unite the thermodynamics and dynam- 
ics of supercooled liquids within one conceptual frame- 
work, and this is the motivation of the present paper. 
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In the past, the two-state model by Angell and Rao^ 
led to configurational entropies that are in quite good 
agreement with experimental data^ The model sug- 
gests that the thermodynamics of super-cooled liquids 
can be described by an ensemble of non-interacting two- 
state excitations^ 5 - each creating an excess of entropy^ 
A recent extension of the models considered two Gaus- 
sian manifolds of levels (2G model) instead of two dis- 
crete states of the original two-state model. This vari- 
ant of the random energy model, conceptually related to 
Basslcr's random-walk picture ) 50 ' 51 allowed an accurate 
fit of the laboratory and simulation data for heat capac- 
ities and configurational entropies. Here, we present a 
somewhat simplified version of this model that consid- 
ers excitations from the single-energy level of the ideal 
glass to a Gaussian manifold of configurationally excited 
states (1G model). The thermodynamic analysis is then 
further used as a basis for a dynamic model of configu- 
rational excitations. 




FIG. 1: Excitable units of a glass-former in real space. 



Our model of dynamics of viscous liquids follows the 
philosophy of the energy models of activation in that 
it considers the probability of accumulating kinetic en- 
ergy sufficient to lift a "structural element" (or "excitable 
unit" , sec Fig. [T]) to an energy level corresponding to the 
activated state of the high-temperature relaxation. The 
kinetic energy supplied by the surroundings of a given 
excitable unit is a fluctuating variable with a fluctua- 
tion width related to the ruggedness of the landscape 
through the configurational heat capacity. The average 
of the fluctuations of the kinetic energy results in a rela- 
tion for the relaxation time which gives an increase of the 
activation energy with decreasing temperature in terms 
of the configurational heat capacity, rather than the con- 
figurational entropy of the Adam-Gibbs theory. The new 
relation gives an account of experimental dielectric relax- 
ation comparable to the Adam-Gibbs formula, but more 
acceptable in several ways to be discussed below. 



II. THERMODYNAMICS OF 
CONFIGURATIONAL EXCITATIONS 

A. Formulation of the model 

We conceive an excitatio n 53 ' 54 ' 57 to be a local increase 
in potential energy that results from a collisional redistri- 
bution of kinetic energy amongst some minimum group 
of "rearrangeable units" of the liquid structure. In this 
redistribution a "unit" (Fig. [T|) that undergoes an ex- 
cessively anharmonic vibrational displacement can get 
trapped by the motions of neighbors that act to prevent 
a return of the unit to its original center of oscillation, 
such that some kinetic energy will be lost from the vibra- 
tional manifold and stored in the configurational mani- 
fold. It is the constant and repeated exchange of energy 
between these manifolds that is the essence of configura- 
tional equilibration. The ability to store potential energy 
by this mechanism determines the configurational heat 
capacity. 

The mathematical realization of this concept assumes 
that a liquid can be divided in real space into excitable 
units that follow the statistics of independent entities. 
Since these units form a continuous dense liquid phase 
and should be interacting, their actual interactions are 
represented by mean-field parameters of excitation en- 
ergy eo, excess entropy So, and excess volume v . Each 
excitable unit is a truly microscopic object as small as 
a fragment of a molecule or a single chemical bond of 
a network glass. The low-energy state is identified with 
the ideal glass, while excitations associated with molec- 
ular motions belong to a Gaussian manifold of energies 
with the width a (Fig. [2]). The energetic disorder arises 
from the local disordering field and packing restrictions 
around a given excitable unit. 

The real-space model of two-state excitations is next 
projected onto Goldstein's configuration space of a liquid 
at constant pressure.— This configuration space can be 
separated into basins of attraction characterized by the 
minimum depth 0. We can write the excess Gibbs energy 
relative to the energy of the ideal-glass state (superscript 
"ex") as a double sum over the basin energies and the 
fraction of excited units x out of N excitable units: 



-g ex JV/T 



-4>N/T 



E < 

0<:c<l 



(7) 




AH 



Here, P(4>, x) is the distribution of minimum energies in 
configuration space obtained by projecting the excitation 
energy x(eo + Pvq + Se) on the Gaussian manifold 



P{4>,x)-- 

where G(Se) is a Gaussian distribution 
G(Se) cxexp [-(6e) 2 /2<t 2 



x(e + Pv Q + 5e)]G(5e)dSe, (9) 

(10) 
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FIG. 2: Configurational excitations and flow/relaxation 
events in the 1G model. The thermodynamics of supercooled 
liquids is represented in real space by excitations from a single 
energy level of the ideal glass into a Gaussian manifold of con- 
figurationally excited states. Each excitation lifts the energy 
of an excitable unit [Fig.fT] by energy eo(x) [Eq. (|22|l ] and en- 
tropy so . The distribution of these states is much broader for 
fragile liquids than for intermediate and strong liquids. Re- 
laxation is described by the dynamics of activated transitions 
in configuration space from basin minima of configurationally 
excited states to a common energy level eo above the top of 
the energy landscape. The activation barrier eo + A<6 is a 
sum of the high-temperature activation energy eo and the 
difference A<f> of the average energy of the basin minima from 
the high-temperature plateau. 



All energies here and below are in K, entropies and heat 
capacities are in units of fee- 

The ratio of vibrational-rotational partition functions 
QQ; e in the ground (superscript "g") and excited (super- 
script "e") states can be absorbed into the excitation 
entropy 



s a = \n[Ql/Ql 



(11) 



which is composed of the harmonic vibrational contribu- 
tion, sjj, and a configurational contribution, Sq. The vi- 
brational excitation entropy is related to the excess den- 
sity of states of low- frequency vibrational modes near and 
below the boson peak^ 



(12) 



In Eq. (|12p , the sum runs over the vibrational frequencies 
uj (eigenvalues of the Hessian matrix) with the densities 
of vibrational states in the ground and excited states g->j e . 

The thermodynamic limit TV — * oo transforms the en- 
tropy s((j),x) in Eq. ||7J) into a sum of the ideal mixing 
entropy, sq(x), and a Gaussian termS 



s(4>,x) = s (x) - 



[</> - x{e + Pv )f 
2x 2 a 2 



(13) 



sq(x) = xsq — x ln(x) — (1 — x) ln(l — x). 



(14) 



The sum over x in Eq. () 1 1 [) is determined by its largest 
summand at x = x{4>). One then arrives at the land- 
scape thermodynamics in which the thermodynamic ob- 
servables are determined by the excess free energy func- 
tion depending on <j> (we omit the dependence on P for 
brevity) 



3 cx (0)=0 + e - h (0)-T S cx (0), 



(15) 



where e| x 



h (cA) is the energy related to the anharmonicity 
effects not included in the harmonic approximation (en- 
tropy from anharmonicity is small as indicated by com- 
puter simulations 3 ^). The excess free energy g cx ((j>) is 
composed of the configurational and vibrational parts 



where 



c (0) =g c {4>) +g v (<t>), 



(16) 



(17) 



One can alternatively consider the sum over <p in Eq. 
([7]) that maximizes at the average basin energy (4>(x)). 
The (partial) Gibbs energy can be considered as a func- 
tion of the population: 



g cx (x) = (^(x)) 



e^ h ((<M*)» 



Ts((cj>(x)),x). (18) 



The minimum of g ox (x) gives the thermodynamic Gibbs 
energy g QX in Eq. The formulation in terms of g ex (x) 
is thermodynamically equivalent to the landscape ther- 
modynamics in terms of g ex (4>) since the thermodynamic 
Gibbs energy g cx is achieved at the largest summand in 
both x and (f> in the double sum in Eq. (J7J). We will, 
however, obtain both g cx (4>) and g cx (x) in order to gain 
better insight into the physics of the model. 

The excitation energy eo can, to the first approxima- 
tion, be considered as independent of temperature. The 
situation is quite different with the Gaussian width a 2 . 
The energy of the localized excited state is randomized by 
interactions with the thermal motions of the liquid which 
arc not quenched and therefore affected by temperature. 
The fluctuation-dissipation theorem then requires that 
a 2 = 2XT scales linearly with temperature ; 62 ' 63 as does 
the mean-square displacement of a classical harmonic os- 
cillator, (x 2 ) cx T . Here, A is the trapping energy or 
the energy of stabilization of an excitation by the dis- 
order of the medium in which it exists. 50 ' 64 Once this 
real space Gaussian width is substituted into Eqs. (JSj)- 
()10p . it results in an approximately linear temperature 
scaling of ((5(f)) 2 ) (which is more complex because of a 
generally non-parabolic form of g cx (4>) and temperature 
dependence of the population x^-) . 

The notion of the linear (cx T) temperature depen- 
dence of ((<5c/)) 2 ) (which we will verify below) is one of the 
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central components of the present model^ The energy 
landscape of the system, i.e. energy as a function of 3N 
coordinates of the molecules making up the liquid, is de- 
termined by intermolecular interactions and is expected 
to be weakly temperature dependent at constant volume 
of the liquid. However, when the manifold of all possible 
states in 37V-space is projected onto one single coordinate 
of the basin energy 0, the distribution of 0, given in terms 
of the (partial) free energy <7 ex (<^), gains temperature de- 
pendence. Not only the first moment of this distribution, 
the average basin energy, is temperature dependent, as 
indeed described by random-energy models^ but essen- 
tially all higher moments arc temperature-dependent as 
well. The random-energy model, originally developed for 
spin glasses,— assumes that the width of the Gaussian 
distribution of random spin configurations is independent 
of temperature. This assumption is well justified for sys- 
tems with quenched disorder, but probably not as well 
for liquids in a metastable (slow nucleation) equilibrium. 
Likewise the original Stillinger- Weber formulations! as- 
sumed that temperature affects only the average energy 
(cf)) in the form of "descending into the landscape" but 
not any higher moments of the distribution. The simula- 
tion evidence on this matter is insufficient and somewhat 
controversial. While some simulations of small ensem- 
bles of binary Lennard- Jones (LJ) fluids give clear indi- 
cation of an approximately linear dependence of ((<5</>) 2 ) 
on T ) 7 i 8 i 68 i 69 simulations of larger systems give virtually 
constant width.— This distinction is not accidental. The 
width of the distribution of inherent structures scales as 
1/y/N and needs to be measured on small ensembles. 

The ideal-glass state should in principle involve ran- 
domness (simulations for network liquids show a finite 
configurational entropy at the cutoff energy 7 -^) and the 
previous version of the model^ 7 - (2G model) assumed 
Gaussian distributions for both the ideal glass energies 
and the excited configurations with the widths of = 
2ksTXi (i = 1,2). However, as we noted in Ref . l57l. the 
fit of the 2G model to experimental heat capacities and 
configurational entropies resulted in a small and almost 
constant Ai ~ 15 — 40 K for the ideal-glass state. The 
application of the model to a more extensive list of glass- 
forming liquids performed in this paper has shown that 
Ai can be set equal to zero without sacrificing the quality 
of the fit. This will be the model adopted here. This ver- 
sion of the model, with only one Gaussian manifold for 
the configurationally excited states, will be referred to as 
the 1G model. The random energy statistics we consider 
here have much in common with the model of protein 
folding proposed by Bryngclson and Wolyne o 59 ! 71 and the 
equations we derive share features with the earlier coop- 
erative two-state model of Strassler and Kitted (which 
is a forerunner of the two species non-ideal liquid model 
of Rapoportf 7 -^ the "two liquids" model of Aptekar— and 
Ponyatovsky and co-workers,— and the cooperative de- 
fects models of GranatoSi and of Angcll-Moynihan 7 - 7 -) . 

Now we turn to the thermodynamics of configurational 
excitations. The excess Gibbs energy <? ex (</>) minimizes 



at the average basin energy 

(0) = x(e Q + Pv ) ~ 2x 2 A(l + deZJd<j>). (19) 

In what follows we will neglect the generally unknown 
derivative <9e° x h /<90. This approximation is expected to 
be accurate at low temperatures close to T g , but will fail 
at higher temperatures above the onset temperature at 
which the system starts to descent into the energy land- 
scape. In this approximation, the excited-state popula- 
tion is defined by the self-consistent equation 



1 



where the free energy per excitable unit is 
9o(x) = eo(z) + Pvo - Ts . 



(20) 



(21) 



Because of the energetic disorder of the exited states the 
actual excitation energy is lowered by twice the trapping 
energy 2A . 57 ' 78 However, only configurationally excited 
states can "solvate" (stabilize) the excitation, and the 
stabilization energy is proportional to the population x of 
the excited states. The effective excitation energy eo(x) 
in Eq. (f2Tj) and Fig. [2] then becomes 



e (x) = eo + Pvq — 2xA, 



(22) 



where the factor of two comes from counting all inter- 
actions of a given unit with the rest of the ensemble of 
configurational excitations. 

By relating x to the spin variable a = 2x — 1 Eq. ([20]) 
can be brought to the form usually considered by models 
of ferromagnetismZ 9 - 



a = tanh 



Act 
2T 



So 



A 



2T 



(23) 



where the excess Gibbs energy of configurational excita- 
tions is 



go = £o + Pvq - Ts . 



(24) 



At so = and A = eo + Pvq Eq. (|23)) transforms into 
the Weiss formula for spontaneous magnetization 7 - 9 - with 
Act/2 playing the role of the effective field of the magnetic 
moments. 

The self-consistent equation for the population of con- 
figurational excitations [Eq. (|20|) ] bears some similarity 
with the results of previous studies minimizing the mean- 
field free-energy functionals of liquid-state theories. The 
density of the liquid p(r) can then be found from the self 



consistent equatio: 



,80,81 



p(r) = gexp ( / c(r — r')p(r')dr' 



(25) 



where q is the activity and c(r) is the direct correlation 
function. Solution of this equation predicts a first-order 
transition to aperiodic crystal.— MD simulations, how- 
ever, indicate^ 2 - that a kinetic transition happens before 
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TABLE I: Number of excitable units ("beads") per mole of a glass-former. The last column represents the effective number of 
classical oscillators in a molecule (mole) of a substance at T g . 



Substance 


Privalko^ 


Takedai^ 


Moynihan^ 


Stevenson^ 


CW(T g )/3 


Inorganics 












Se 


1 


1 


1 


0.9 


1 


ZnCl 2 


1 






1.2 


3 


Aromatics 












Toluene 


4 




1 


2.7 


3.2 


o-terphenyl (OTP) 


9 




2 


3.7 




Paraffinics 












2-methylpentane 


6 




3 


3.8 


3.2 


Alcohols 












methanol 


2 


2 




1.3 


1.76 


glycerol 


6 


6 


7 


4.5 


3.5 


Hydrates 












Ca(N0 3 ) 2 -4H 2 


13 






7.0 


11.7 



the thermodynamic transition is reached, in qualitative 
agreement with our analysis of experimental data (see 
below) . 

The combination of the average energy [Eq. (flTj)) ] and 
excess entropy [Eqs. (|T^|) and (fT5)) ] as functions of popu- 
lation yields g cx (x) in Eq. (fT8j) 

+ T[xlnx+ (1 -x)ln(l - x)] . ' 

Except for the anharmonic correction, this equation has 
been derived in many previous publicat ions i 72 ' 73 > 74 ' 75 1 83 
Minimization of g ex (x) with respect to x (neglecting the 
derivative de^ h /dx) leads to Eq. (|20|) . The parameter A 
then plays the role of the average energy of interaction 
between the configurational excitations. Consequently, 
the quadratic in x term in Eqs. (|19[) and (|26p . originat- 
ing from the energy randomness in our model, is equiva- 
lent to direct mean-field (Bragg- Williams 7 ^) interaction 
between the excited units. In other words, randomness 
is effectively equivalent to attraction when separate units 
independently seek the same satisfactory configuration^ 
Notice that derivation of Eq. (j2l))) requires the explicit ac- 
count of the linear temperature scaling of the real-space 
Gaussian width a 2 in Eq. (fTU|) . The assumption of a 
temperature-independent width would result in a 1/T 
scaling of the energy term in <? ox (x) quadratic in x. 

Since the full calculation of the excess thermodynamics 
of the supercooled liquid over its crystal is too a complex 
task even for phenomenological models we next assume 



that the excess entropy of configurational excitations over 
the ideal-glass state, s ex (T), gives the excess entropy of 
the liquid over its crystal. We will use this excess en- 
tropy for the rest of our thermodynamic analysis warning 
at this point that this entropy is not thermodynamically 
consistent with the excess Gibbs energies in Eqs. (fl5]) 
and dM}, i.e. s 0X (T) ± - (dg cx (T)/dT) p . A thermody- 
namically consistent g ex (T) can of course be obtained by 
temperature integration of s ex (T). 

According to our derivation above, s ex (T) is the sum 
of the configurational entropy s c (T) and the vibrational 
entropy x(T)sq: 

s cx (T) = s(((f>(x(T))) , x(T)) = s c (x(T),T)+x(T)sl 

(27) 

where (0(x(T))) and x(T) are given by Eqs. (fTTJ)) and 
PI!]) , respectively. The configurational entropy thus ac- 
commodates only that part of sq which is not related to 
the change in the vibrational density of states. It is given 
as a sum of the ideal mixture entropy Sq(x) and the fluc- 
tuation entropy — x 2 \/T consequent on the shrinking of 
the Gaussian width with decreasing temperature'^ 7 - 

s c (x,T) = s c (x)~x 2 \/T, (28) 

where 

Sq(x) = xSq — xlnx — (1 — x) ln(l - x). (29) 



B. Application to experimental data 



We now proceed to applying the model to laboratory 
and simulation data for supercooled molecular liquids. 
In the former case, excess entropy As of the liquid over 



the crystal is identified with the excess entropy over the 
ideal-glass state s ex (T) [Eq. ([27]) ]. In the case of sim- 
ulations, conjugate gradient minimization of simulated 
trajectories allows sampling of inherent structures, and 
equations for the configurational component of the ex- 
cess entropy will be used [Eqs. (f28|) and (|29|) ]. The two 
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TABLE II: Best-fit parameters of the 1G model to experimental excess entropies and heat capacities. Also shown are the 
experimental Kauzmann temperature T^ p obtained by extrapolating experimental entropies to zero (Ref. IHih and by using Eq. 
(|38[) . experimental VFT temperature To, and the thermodynamic Kauzmann temperature Tk calculated from Eq. Q using 
the 1G excess entropy from Eqs. (|27p - (|29|l . The temperature of liquid-liquid transition Tll is calculated from Eq. (|34[1 . T g is 
the experimental glass transition temperature. All energies and temperatures are in K, the entropy so is in fce units. 



Substance 


a, b 

m ' 


z c 


eo 


A 


so 


rpCXp 

K 


To" 


rpCXp d 

1 K 


Tll 


T^ 


T g 


Fragile liquids 


Toluene 


W5 } 


1 


2171 


1020 


10.2 


100 


96.5 C 


100 


113 


108 


117 


D,L-propene carbonate (PC) 


104 


1 


2921 


1383 


10.7 


129 




129 


144 


123 


156 


o-terphenyl (OTP) 


81 


2 


3576 


1686 


8.3 


204 


202.4 


206 


228 


170 


246 


2-methyltetrahydrofurane (MTHF) 


65 


2 


899 


414 


5.9 


69 


70 


70 


82 


50 


91 


Salol 


63 


3 


2070 


988 


5.6 


175 


175 


176 


193 


115 


220 












Intermediate 


liquids 










3-bromopentane (3BP) 


53 


4 


348 


20 


2.0 


84 


83 






75 


108 


Glycerol 


53 


7 


738 


5.0 


1.6 


137 


130 9 






59 


190 


n-propanol (nPOH) 


35 


2 


406 


22 


2.8 


72 


70 






30 


96 



"Steepness fragility index, Eq. J33P . 

6 Taken from Ref. [27] unless indicated otherwise. 

c From Ref. |H 

Calculated from Eq. 1551 . 

c Obtained as the temperature of crossing zero for the excess en- 
tropy calculated in the 1G model. 
^Ref. 
»Rcf. 



sets of equations are formally equivalent when the overall 
excitation entropy sq is used for the excess data and only 
its configurational component sjj is used for the statistics 
of inherent structures. 

No progress can be made without recognizing that 
the molar quantity, heat capacity, determined in a lab- 
oratory measurement contains contributions from z "re- 
arrangeable subunits" of the molecule which are the mi- 
croscopic dynamic elements of the system. This is most 
obvious in the case of chain molecule systems like sele- 
nium or like 4-methyl nonane^ where each methyl group 
constitutes one "bead" according to Wundcrlich's origi- 
nal description.— The number of subunits or beads per 
molecule is not always evident. In selenium the atom 
is clearly the "bead" but in the glass-former 9-bromo 
phenanthrene^l the whole C20 "raft" is rigid and can 
only rearrange as a single unit. 

Stevenson and Wolynes have suggested an operational 
approach to determining z, based on the entropy of fu- 
sion relative to that of the simple Lcnnard Jones (LJ) 
system^ While this method does not take account of the 
fact that the LJ fusion entropy is determined at a temper- 
ature where the liquid is enormously more fluid than the 
glass-formers at their melting points (entropy per bead 
1.68 ks compared to Wundcrlich's 1.36 kp, which leads 
to gross overestimates in the case of low fusion entropy 
network liquids, e.g. ZnC^), it does give a measure that 



is roughly consistent with others. These range from the 
original Wundcrlicli 86 and Privalkc4£ estimates, through 
Takeda et al. based on structural arguments^ down to 
Moynihan and Angcll whose estimates were based on best 
fitting of the excess entropy to an excitations models 
Some comparisons are provided in Table HI 

The separation of a molecule into z statistically inde- 
pendent subunits neglects the finite correlation length of 
the disordering field. Indeed, if two units are within the 
field's correlation length, they contribute to the statistics 
of basins as a single unit. The number of units z is thus 
an effective parameter necessarily smaller or equal to the 
number of conformationally distinct units (cf. the differ- 
ence between Privalko's conformational and Stevenson's 
effective numbers in Table HI)) . Because of its effective na- 
ture, z can be represented by fractional numbers as was 
done by Stevenson and Wolynes (Table UJ. Although this 
approach provides more flexibility in fitting the experi- 
ment, we will follow here Takeda et al£l and Moynihan 
and Angell^ and use integral numbers for z. In order 
to distinguish between thermodynamic properties refer- 
ring to excitable units and molecules (or generally moles), 
we will use lower-case letters for the former and upper- 
case letters for the latter. Lower-case and upper-case 
extensive variables are connected through z, e.g. for the 
constant-pressure heat capacity, Cp = zcp. 



We will apply the model to constant-pressure data and 
consider the excitation energy eo as an adjustable param- 



eter. Therefore, for the rest of our analysis, the term 
Pvq in the excess Gibbs energy [Eq. (j24|) ] is fused into 
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erj. Therefore, given the number of excitable units per 
molecule z has been specified, the 1G model contains 
three model parameters, e , A, and s . We have tested 
the model for a set of glass-formers identified as fragile 
(strong super- Arrhcnius kinetics, from toluene to salol) 
and intermediate (weak super- Arrhenius kinetics, from 
3-bromopentane to n-propanol) liquids. Results shown 
in Fig. [3] and listed in Table HT1 were obtained from a si- 
multaneous fit of the model to excess heat capacities and 
entropies from Ref . l54l. bound by the constraint e (x) > 
[Eq. (|22p] required for the mechanical stability of the 
ideal-glass state. 

The experimental excess entropies are calculated from 
the constant-pressure heat capacity Cp K (T) and the fu- 
sion entropy ASf us according to the relation 



S ex (T) = A5 fua + 



Cp K (T')(dT'/T'), 



(30) 



T fua 



where Tf us is the fusion temperature and 5 ex (T) = 
zs cx (T). The number of excitable units per molecule 
(mole) z was taken from Moynihan and Angell^ and 
varied additionally to find the best-fit integral numbers. 
Equal quality fits can be obtained in some cases with 
different numbers z, e.g. for toluene z = 1 and z = 2 
can be adopted. In that latter case, z = 1 was taken to 
maintain the consistency of parameter values with other 
fragile liquids (Table [TTJ) . The choice of z here does not 
affect our qualitative conclusions discussed below. 



shift (which is an analog of the trapping energy for elec- 
tronic transitions) on temperature^ With this assump- 
tion one gets: 



x 2 X 
~T~ 



x(l 



e(x)(e(x) - 2xX) 
*T 2 (l-x(l-x)(2X/T))' 



(31) 



At A = this equation reduces to Schottky's heat ca- 
pacity c c p = x(l — x)(eo/T) 2 , which is further reduced 
to the Hirai-Eyring equation* 9 ^ proposed on the basis of 
transition-state ideas, in the limit x <C 1. For fragile 
liquids, x ~ 1, as we show below, and the heat capacity 
becomes: 



cf? = X/T. 



(32) 



Notice that Eq. (j31"j) anticipates a Curie- type, (T— T c ) , 
divergence of the heat capacity at the critical tempera- 
ture defined by the equation T c = 2Xx(T c )(\ — x(T c )). 
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FIG. 3: Excess entropy (a) and excess heat capacity (b) 
for some of the supercooled liquids listed in Table [TT] (per 
molecule, in units of kn). The thin lines refer to fits to the 
1G model, the thick lines refer to the experimental data.— 

The constant-pressure heat capacity, used to fit the 
experimental data, was obtained by the direct differ- 
entiation of the excess entropy constrained by the as- 
sumption that the trapping energy A is independent of 
temperature. This assumption is supported by spectro- 
scopic studies showing weak dependence of the Stokes 



FIG. 4: Excess entropy (solid lines) and its ideal-mixture com- 
ponent zsq(x) (dashed lines) [Eq. (|14|l ] vs temperature. The 
dashed and solid lines coincide on the scale of the plot for 
glycerol. 

The fragility of a glass-former is often characterized by 
the steepness index^ 

m. = d\ogr/d(T g /T)\ T=Tg (33) 

(listed in Table [TTJ) or by its thermodynamic 
equivalent! 92 ! 93 What appears to be a smooth transi- 
tion from fragile salol to intermediate 3-brombenzene 
according to the steepness index in fact corresponds to 
a drastic decrease in the trapping energy by a factor of 
about 50. Low values of A turn out to be characteristic 
of all intermediate liquids in Table [TTJ The result is a 
profound change in the relative importance of the ideal 
mixing and Gaussian terms in the excess entropy. 

For fragile liquids, the ideal mixing entropy [Eq. ([T4"]) 
or Eqs. (f2"T|) and (f2U)) ] is large and is almost constant 
in the whole temperature range of a supercooled liquid 
T g < T < T m (toluene in Fig. gj). This is a reflection 
of the fact that the population of the high-energy state, 
driven by the excitation entropy so, is close to unity 9 ^ 
and almost independent of temperature. The excess ideal 
mixing entropy is finally lost in a first-order transition to 
a low-entropy liquid state at a temperature Tll which, 
according to fits to experimental data, lies between the 
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FIG. 5: Temperature dependence of the average basin energy 
(a) and excess entropy (b) of fragile glass-formers. The al- 
most flat dependence of the average energy on temperature 
is terminated by discontinuous first-order phase transition at 
temperature Tll- At this temperature the excess entropy, as 
predicted by the 1G model [solid line in (b)], also shows a 
discontinuous drop to a nonzero value which decays to zero 
at the thermodynamic Kauzmann temperature Tk- The bro- 
ken line in (b) illustrates the experimental interpolation of 
the high-temperature 1/T law ending up at the experimental 
Kauzmann temperature T^ p . The excess entropy is equal to 
the fusion entropy A5f us at the melting temperature T m . The 
dashed lines in (a) indicate the narrowing of the Gaussian dis- 
tribution of excitation energies with decreasing temperature. 

Kauzmann temperature Tk and T g (see Table [TTJ) . The 
temperature of the equilibrium phase transition is defined 
in terms of the model parameters as 

T LL = (34) 
so 

This first-order transition (Fig. EH) gives rise to a peak 
in the heat capacity which is cut off by the kinetic glass 
transition at T g . The observed drop of the heat capacity 
is related to the loss of ergodicity which we do not con- 
sider here. The entropy does not drop to zero at Tll but 
in all cases becomes small, crossing zero at the Kauzmann 
temperature Tk much below the experimental tempera- 
ture T^ p (Table [TTJ) . A second order phase transition 
can be realized at the critical point when the excitation 
energy gap is related to the trapping energy by the fol- 
lowing equation^ 

e = A + As /2. (35) 

However, when this restriction is imposed on the param- 
eters, the model fails to fit the experimental data. 



The physics is quite different for intermediate liquids 
for which the ideal mixture term provides the main part 
of s 0X (T). The excess entropy then decreases with tem- 
perature due to decreasing population of the excited 
state, which cannot stay at a high value because of a 
smaller entropy gain sq compared to fragile liquids (glyc- 
erol in Fig. 3]). The entropy smoothly decreases without a 
discontinuity related to thermodynamic phase transition 
(see Fig. [5]3 for the illustration of various temperatures 
used in the present thermodynamic analysis). 

The disappearance of the first-order transition for non- 
fragile (strong and intermediate) liquids is related to two 
critical parameters, critical temperature T c and critical 
excitation entropy, sq c - A first-order phase transition is 
possible for temperatures below T c and entropies higher 
than sq c £Z- 

T < T c , s > s 0c , (36) 

where 

T c = A/2, s 0c = 2e /A - 2. (37) 

At least one of two inequalities in Eq. (|36|) is violated for 
non-fragile liquids as a result of low excitation entropies 
So and trapping energies A. 

In parallel to the excess entropy, the average basin en- 
ergy (0) shows qualitatively different temperature behav- 
ior for fragile and intermediate liquids. With tempera- 
ture decreasing, ((f)) from Eq. (fTD)) starts to dip as 1/T 
from a high-temperature plateau and then inflects into 
an exponential temperature dependence oc exp(— eo/T) 
for e(x)/T » s . Both the 1/T decay at relatively high 
temperature a 30 ' 33 i 68 ' 95 and exponential decay at low tem- 
peratures (for a model network fluid2£) have been ob- 
served in simulations. In our model, this pattern de- 
scribes non-fragile liquids. For fragile liquids, the excited 
state population is almost constant, x ~ 1, in the entire 
range of experimentally accessible temperatures down to 
the glass transition. The average basin energy is con- 
stant as well, ((f)) = eo — 2A. For infinitely slow cooling 
(and rare cases like triphenyl phosphite (TPP)), the pop- 
ulation sharply changes at the liquid-liquid transition re- 
sulting in a discontinuous dip of the average energy ((f)) 
from its plateau value. 

The temperature dependence of the excess entropy of 
fragile liquids above T g is mostly determined by the 1/T 
decay of the second, fluctuation term in Eq. ([25]). The 
result is the overall temperature dependence in the form 
of Eq. ^ with 5*o = zsq. The experimental Kauzmann 
temperature is then obtained by extrapolating Eq. © to 
zero entropy, which leads, in terms of the model param- 
eters, to the following relation 

T£ xp = -. (38) 
so 

Equation (j3"5|) holds quite well for the fit parameters in 
Table HIl In addition, the constant-pressure heat capacity 
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FIG. 6: Configurational entropy (a), configurational constant- 
volume heat capacity (b), and the fitted excitation parameters 
(c,d) for LW o-terphenyl (z = 2 is used as in the case of 
laboratory o-terphenyl, Table HI)) . The points in (a) and (b) 
are the simulation data from Ref.[3(] at the densities indicated 
in (c). The dotted lines in (c,d) connect the points. In (c), A 
is indistinguishable from eo /2 on the scale of the plot for two 
lowest densities. The lines in (a) and (b) are simultaneous fits 
of the simulated configurational entropy and heat capacity to 
the 1G model; the dashed and dash-dotted lines are used for 
two lowest densities with high fragility to distinguish them 
from higher densities with intermediate fragility. 



scales as 1/T [Eq. (52J) ]. Therefore, parameter A can be 
set by the heat capacity at the glass transition, 



A 



T g c<?{T g ). 



(39) 



The fit of experimental data for fragile liquids shows 
that eo is close to 2A and in fact can be put equal to 2A 
without sacrificing the accuracy of the fit. As a result, 
the thermodynamics of fragile liquids is defined by the 
following relations often used cmpiricaii y 21 i 27 ' 84 i 85 ' 97 



(T) = c5?(T„) (T g /T^ p Tg/T) , 



TO(T S /T). 



(40) 



We note that random energy models, which do not antic- 
ipate temperature variation of the width of basin energy 
distribution^ result in T~ 2 scaling of the heat capacity 
inconsistent with Eq. (|40|) . 

It is obviously significant for our discussion to com- 
pare the results of fitting the model to laboratory ex- 
cess entropies with analogous fits to configurational en- 
tropies from simulations. The Lewis and Wahnstrom 
(LW) model of o-terphenyl studied in Rcfs . [30l and [4ll pro- 
vides us with such an opportunity. Figure [J5] shows the fit 
of the 1G model to the configurational entropies and heat 
capacities of LW o-tcrphenyl at five different densities.— 
The number z = 2 is maintained equal to the analysis 
of experimental data to access the fraction of configura- 
tional excitation entropy in sq. Two main results follow 
from the fit: (i) the 1G model predicts the existence of 
a critical A-singularity below T g for two lowest densities 
and its disappearance at higher densities, (ii) the config- 
urational component of the excitation entropy following 
from the fit, Sq ~ 2, is significantly lower than sq ~ 8 
obtained for laboratory o-terphenyl ( Table HT)l . The ap- 
pearance of the divergence at lower densities is the result 



of very close fulfillment of Eq. (|3"5|) for the fitting pa- 
rameters, which is remarkable, given that the excitation 
parameters are freely varied in the fit. The magnitude of 
the fitted excitation entropy suggests that a major por- 
tion of so arises from a change of the vibrational density 
of stateS ) 37 ' 38 i 77 as we have also learned from comparing 
the laboratory and simulation data for water (see below). 
Our result is close to Goldstein's estimate of 28% for the 
fraction of configurational component in the excess en- 
tropy of o-terphenyl.— Goldstein's calculation was based 
on the comparison of excess entropies of quenched and 
annealed glasses at T g , ASQA{T g ), and at K, ASqa(0)- 
The ratio ASqa{0)/ ASQA(T g ), which is equal to sg/s 
in the 1G model, gives the configurational fraction of the 
excess entropy. 

The fit of 1G thermodynamic model to experimental 
data (Table [TT| shows that the thermodynamic Kauz- 
mann temperature Tk [Eq. ^] is significantly lower than 
the experimental Kauzmann temperature T^ p [Eq. ([5])]. 
This implies that the relaxation time in the Adam-Gibbs 
relation does not diverge at T™ p ~ T and the link be- 
tween the thermodynamics and dynamics might be more 
complex. The dynamic extension of the 1G model pre- 
sented below places the emphasis on the configurational 
heat capacity, instead of the configurational entropy, as 
the main reason for super- Arrhenius dynamics in fragile 
liquids. 



III. DYNAMICS OF CONFIGURATIONAL 
EXCITATIONS: ENTHALPIC DRIVING FORCE 

A. Formulation of the model 

Here we describe a dynamic model extending the ther- 
modynamic analysis of Sec. |TTJ Our development starts 
with the assumption common to all energy-trap^^ and 
cncrgy-diffusio n 50 i 51 i 52 models that non-reversible events 
of viscous flow and diffusion^ occur by exciting some 
states within the liquid to a common energy level Eq, 
which is higher than the top of the energy landscap o 61 i 98 
(Fig. [2]). These excitations occur by absorbing kinetic 
energy by an excitable unit (bead) from the surrounding 
liquid. 

We will next assume that only unjammed configura- 
tionally excited units will participate in activated events. 
Even when a sufficient amount of energy has been ac- 
cumulated at a given unit, relaxation event may require 
facilitation from other units. This is particularly clear in 
a case of a molecule composed of z units (beads). One 
could imagine that e.g. translational relaxation of such 
molecule would require excitation of all z units, although 
relaxation of conformationally flexible molecules can also 
proceed in a diffusive way, as a sequence of low-amplitude 
motions of consecutive units. Because of the assumed low 
amplitude of the motions involved, we will not distinguish 
between pairs of units within the molecule and pairs of 
units belonging to different molecules. 
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hot region 




configuration 



excitation 

FIG. 7: Hot region of N c configurationally excited units 
(hatched circles) with kinetic energy accumulated by transfer 
of m "quanta" of the kinetic energy ex from the surrounding 
liquid. 



Each low-amplitude step is terminated by a transfor- 
mation of the accumulated kinetic energy into a small 
structural change with a higher potential energy and thus 
higher fictive temperature. A single relaxation event re- 
quires a sequence of such locking steps by the end of 
which there appears a "hot" (in fictive temperature) re- 
gion within the liquid which can then subsequently relax 
to a new configuration (Fig. [7]). In this scheme the over- 
all dynamics become hierarchical in character, the initial 
step being the most probable and hence the shortest in 
time. This picture is an integration of energy spikes seen 
in simulations of Heuer and co- worker s 95 ! 98 with dynam- 
ical hierarchy of Palmer et al£2- (see below) . 

The kinetic energy becomes increasingly scarce at low 
temperatures. The creation of a hot region will occur 
by pulling the kinetic energy from a growing number 
of neighboring molecules, leading to the creation of a 
hot island in a sea of kinctically frozen molecules (on 
the time scale of heat transport). This picture bears 
some similarity to the entropy-rich droplet enveloped by 
the entropy-frozen environment described by Lubchenko 
and Wolynes^ and by Bouchaud and Biroli^ How- 
ever, the notion of relaxation proceeding by occurrence 
of hot regions does not anticipate static, thermodynami- 
cally stable structures of the mosaic picture^ and in- 
stead corresponds to the idea of dynamic heterogene- 
ity i.e. the existence of regions of markedly different 
mobilityi 101 i 102 Computer simulations generally support 
this viewi 31 i 103 i 104 Within this general umbrella, dynamic 
heterogeneity can be treated in two distinct ways: as a 
static distribution of relaxation times in regions of vary- 



ing mobilit y 21 i 102 i 105 or in a dynamic fashion as facili- 
tated kinetics of transfer of excitations from mobile par- 
ticles to their neighborsi 5 - 5 - Our formulation will follow 
this latter pathway using the kinetic scheme of hierarchi- 
cal relaxation events advanced by Palmer et al^ 

According to the hierarchically facilitated dynamics, 
step n + 1 happens only when a configuration at step n is 
reached to facilitate the next move^ This idea has been 
instrumental in establishing the conceptual basis for ki- 
netically constrained models^ and, physically, leads to 
dynamical heterogeneity when dynamically cooperative 
regions in the glass are created by a sequence of con- 
strained motions of fastest molecules in the ensemble^ 
In our model, step n is reached when n excitable units 
within a hot island have "blinked" into the excited state 
with the probability 

Pn = (Vti)" = exp[H, (41) 

where t± is the average waiting time for a single unit and 
To is the same as in Eq. $Z§. As a result of a sequence 
of correlated steps, each resulting in excitation of n units 
out of N c units in the hot island, the waiting time of level 
n becomes^ 9 - 

t„ = T exp[^7V„] (42) 

with N n = n(n — l)/2. We can now follow Brey and 
PradosiSZ. to obtain the normalized relaxation function: 

4>{t) = C 1 [E^te^/ro) - E^t/ro)] , (43) 

where E\(z) is the exponential integral function, 

C = M^max, (44) 

and 

N max = N C (N C - l)/2. (45) 




ln(t/x ) 

FIG. 8: (j>(t) from Eq. (|43)l (solid lines) and the closest cor- 
responding KWW relaxation function [dashed lines, Eq. JT}] . 
The values of £ used to plot (j>(t) are indicated in the plot. 

For most practical purposes the relaxation function in 
Eq. (g5J) is indistinguishable (Fig. [8]) from the KWW 
function in Eq. {!]). In the intermediate range of times, 
1 <C t/ro <C e^, which can be very broad since £ ^ 1 for 
real systems, 4>(t) follows the logarithmic decay, cf>(t) ~ 
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FIG. 9: Activation barrier from the average energy (e) to the 
energy level eo above the top of the landscape and fluctua- 
tions of the kinetic energy e_R- in a rugged landscape of fragile 
liquids. The activation barrier eo refers to relaxation in a 
high-temperature liquid, A(fi = (<j>(T)) — (<jf>(oo)} is the change 
of the average minimum energy from the high-temperature 
plateau. The real space regions on the right illustrate fluctu- 
ations in the kinetic energy of the molecules due to variations 
of the basin depths. 

1 — ln(t/ro), common for biopolymersJ^ The average 
relaxation time r = / °° t<f>(t)dt from Eq. (|4"5)) is 

t = (to/C) (e C - 1) - r e c . (46) 

Notice that the KWW function is more flexible than Eq. 
([4*3")) because it involves two free parameters, r/r and (3, 
in contrast to the single parameter £ in Eq. (|4"3")) . More 
complex facilitation rules than the ones used here will 
provide additional parameters and a possibility of real- 
izing the KWW relaxation function^ which has the ad- 
vantage of allowing an approximate correlation between 
the stretch exponent and fragility,!^ 

Each elementary step within the hierarchical sequence 
requires overcoming the activation barrier between the 
average energy of the basin minimum (0) and the com- 
mon energy level E = zeo (Fig. [9]): 

E D (x)=E D -zA<f>(x), (47) 

where Ed is the activation barrier per molecule (mole) 
associated with the activated relaxation in the high- 
temperature liquid^ and 

A0(x) = {x- l)e - 2(x 2 - 1)A (48) 

is the drop of the minimum energy from the high- 
temperature plateau below the onset temperature. 

One can next assume that the kinetic energy neces- 
sary for activation is distributed in "quanta" of thermal 
kinetic energy Ek = zen throughout the liquid, where 
Ek and en refer to a molecule (mole) and excitable 
unit, respectively. This physical picture seems appro- 
priate for describing activated events in disordered ma- 
terials since, according to theories of heat conductivity, 
the quasi-lattice vibrations in glasses are more appro- 
priately described (in the temperature range above ca. 
30 K) as quasi-localized vibrations rather than wave-like 
motionsiiii The momentum exchange (heat transport) 
occurs by diffusional transport of vibrational energy by 



these quasi-localized modes in the vicinity of the boson 
peak i 112 i 113 Activation of one unit then requires accu- 
mulating m = Ed{x)/Ek "quanta" of kinetic energy, 
provided that that unit is not jammed being in the con- 
figurationally excited state. Therefore, the probability of 
activating one unit is equal to the probability of absorb- 
ing m quanta of kinetic energy out of a manifold of xN 
configurational excitations uniformly distributed over N 
units (beads) in the liquid. This type of problem is con- 
sidered in the theory of unimolecular bond dissociation. 
The solution by Kasselii^ gives the probability of com- 
bining at least m quanta of energy at one bond: 

{N + xN -m-l)\(xN)\ 
- m (xN - m)\(N + xN - m)V 1 ' 

For condensed-phase problems one takes the thermo- 
dynamic limit in the above equation, N — > oo, with the 
result 

P> m = (1 + l/x)- n = cxp (-E(x)/E K ) , (50) 
where 

E(x)=E D (x)]n{l + l/x), (51) 

and E D (x) is given by Eq. (|4"7) . 

Because of similar combinatorial rules, not surpris- 
ingly, Eq. (|50|) is analogous to the equation for the prob- 
ability of finding a hole with the volume exceeding some 
critical volume v* 

P(v*)=exp{-jv*/v f ), (52) 

where Vf is the free volume and 7 is a numerical coef- 
ficient. This equation is the key result of free- volume 
models of diffusion and relaxation in glass- formcrs i 17 > 18 
The ideal glass transition is predicted to occur because 
the liquid is supposed to run out of free volume at a fi- 
nite temperature. In contrast to that, the kinetic energy 
in Eq. (|50[) is a fluctuating variable which can approach 
zero for some basins in the distribution sampled at tem- 
perature T (e^ in Fig. [5]), but whose average value is 
proportional to T. 

The combinatorial arguments of the Kassel model envi- 
sion the system as a microcanonical ensemble character- 
ized by the average energy (e) uniformly distributed over 
the sample (Fig. [9}. In order to apply these combinatorial 
rules to a macroscopic liquid, we need to use the micro- 
canonical ensemble characterized by the average energy 
(e) [analogously to the use of microcanonical ensemble to 
calculate the excess entropy in Eqs. ([8]) and (fl~3|) ]. Since 
energy is an extensive variable, this description is equiv- 
alent in the thermodynamic limit to the canonical one in 
the sense that the fluctuations of the total energy can be 
neglected. 

Each unit undergoing excitation to the common level 
eo by collecting kinetic energy from the surrounding 
molecules will find itself in a local disordering field char- 
acteristic of a particular basin of the rugged energy land- 
scape (Fig. [9]) . The kinetic energy available to the sur- 
rounding molecules to excite a given molecule will be 
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a fluctuating variable producing disorder of ex which 
is quenched on the time-scale of momentum relaxation. 
The relaxation time of a single molecule t\ then needs 
to be averaged over possible realizations of the kinetic 
energy 



ti = r e M = t J exp[E(x)/(ze K )]P(e K )de K , 



(53) 



where P{ex) is the distribution of the kinetic energy. 

For the energy landscape characteristic of strong and 
intermediate liquids the canonical narrow distribution of 
basin energies projects itself, in the microcanonical en- 
semble, into a narrow distribution of kinetic energy ej(- 
The distribution function P(ex) in Eq. ([5"3")) can be re- 
placed by a delta function. Assuming that the average 
kinetic energy is equal to Ek = (3/2)T (translations for 
diffusion and viscous flow and rotations for dielectric re- 
laxation), and using Eq. ([5T|) for E{x) in Eq. (|53|) one 
gets for the average relaxation time 

ln(r/r ) = (2/3T)(E D - zA0(x)) In [2 + e -«o+ 6 (*)/r 

(54) 

The relaxation is Arrhcnius at high temperatures 
(e(x)/T <C so an d zAtfi <C Ed)- Two things happen 
when the temperature is lowered. First, the energy gap 
between the average minimum energy (<p m ) and the en- 
ergy level eo starts to increase^ as a + b/T since Acf>(x) 
scales as 1/T right below the onset temperature^ Sec- 
ond, at e(x)/T ^> so the logarithmic term in Eq. ([M]) 
generates a 1/T factor. Overall, the temperature law at 
these low temperatures becomes 



Ht/t ) = E 1 /T 2 + E 2 /T 3 , 



(55) 



which is a linear combination of the Bassler— and 
Litovitaii£ temperature laws. 
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FIG. 10: Test of Eq. ([54} using dielectric relaxation data for 
intermediate liquids listed in Table [TT1 The high-temperature 
activation energy Ed in Eq. (f54)l is equal to 18.4T 9 (glyc- 
erol), 22.6T 9 (nPOH), and 16.0T 9 (3BP)iiiS E D for 3BP was 
obtained from it s em pirical connection to the steepness index 
m given in Ref. 1 1161 . The excited-state population x(T) and 
A</>(T) in Eq. (|54[) are calculated from the 1G model with the 
liquid parameters listed in Table [TTJ The slopes of the linear 
regressions are 1.0 (glycerol), 1.1 (nPOH), and 1.0 (3BP), and 
the intercept gives the physical value of 10 -1,3 s. In case of 
glycerol, the average kinetic energy of 3T has been adopted. 



This sort of iron- Arrhcnius kinetics is capable of de- 
scribing the relaxation of intermediate liquids studied in 
Sec. |TTJ Figure [TO] shows experimental dielectric relax- 
ation tiniest vs the rhs of Eq. ([5"4"|) calculated from the 
1G model. The parameters affecting population x are 
taken from our thermodynamic analysis summarized in 
Table [TTJ Experimental high -temperature activation en- 
ergies Ejj are from Ref. Ill6i The analysis yields straight 
line with the slope 1.1 for n-propanol and a less clear 
linear trend for 3BP with the linear regression slope of 
1.0. For glycerol, the slope of 1.0 is obtained by assum- 
ing the average kinetic energy equal to 3T, which might 
reflect the participation of both rotational and transla- 
tional degrees of freedom in the relaxation process. We 
need to note that vibrational heat capacity of many glass- 
forming liquids does not reach the Dulong and Petit limit 
in its low-temperature portion. Therefore, the average ki- 
netic energy can fall below (3/2)T also becoming a non- 
linear function of temperature. Accounting for this effect, 
which is not considered here, will make the temperature 
dependence in Eq. (|54p even more complex. 

Equation (|54[) is insufficient to account for super- 
Arrhcnius behavior of fragile liquids from Table [TTJ This 
is because fluctuations of the kinetic energy now need to 
be taken into account [Eq. lj5"3"]) ]. Ek can be connected 
to the fluctuation of the energy of inherent structures by 
assuming that the average energy is (E) = (3 /2)T+z ((/>}. 
The instantaneous kinetic energy is then 

E K = (3/2)T - z8<j>, (56) 

where 6(p — <f> — (cf>). The distribution of 8<p is Gaussian, 

W) 2 



P{S(t>) = (2irc c P T 2 ) 1/2 exp 



2c c p T 2 



(57) 



This follows from the bilinear expansion of the enumer- 
ation function s c (4>) m 6<p and the use of two thermody- 
namic identities^ 



(ds c (<t>)/d<f>) NiP = T- 
(c>V(0)/d^ 



>N,P 



(cpT 2 ) 



(58) 



where dp is the configurational heat capacity of the ex- 
citable unit [in contrast to the excess heat capacity in Eq. 
([3"T]) ]. Notice regarding Eq. (|5T[) that Boltzmann statis- 
tics of the basin fluctuations follow from the hyperbolic, 
dp oa 1/T scaling of the configurational heat capacity, 
while 1/T 2 scaling of the random energy models^ leads 
to the temperature-independent statistics. 

The average waiting time of one-unit excitation then 
incorporates the thermodynamic quantity, configura- 
tional heat capacity, into the kinetics through P(S(j>) in 
the integral: 

r 3T/2z 



n/ro 



exp 



E{x) 



3T/2 - z5<j) 



P{8(j})d8(j}. (59) 



A simple estimate of the integral in Eq. can be ob- 
tained by linearly expanding the exponent in 6<j), inte- 
grating over Scj), and reverting to the fractional form. 
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TABLE III: Best- fit parameters of Eq. (|61j) to experimental dielectric relaxation data with To, D, and T' considered as fitting 
parameters. The experimental (superscript "exp") and calculate (superscript "calc") pressure variation of the glass transition 
temperature dT g /dP is given in K/MPa, all temperatures are in K. 



Substance 


log (to /s) 


D 


r 


rj-if a 


rp/b 


dT 9 /dP cxp 


dTg/dP^ 


Toluene 


-15.6 


145 


9.3 


9.8 


8.7 






o-terphenyl (OTP) c 


-14.0 


161 


13.8 


12.2 


13.7 


0.26 d 


0.22 


2-methyltetrahydrofurane (MTHF) 


-13.5 


95.6 


7.9 


5.9 


8.9 






Salol 


-14.4 


204 


11.9 


10.5 


13.3 


0.20 e 


0.26 


3-bromopentane (3BP) 


-13.6 


209 


7.1 




9.4 


0.03 / 




Glycerol 


-14.2 


159 


13.2 




10.0 


0.05 9 


n-propanol (nPOH) 


-13.2 


218 


7.9 




6.5 


0.07 h 


0.05* 



"Calculated for fragile liquids from Eq. I|65|l by applying the ther- 
modynamic fitting parameters from Table [TT1 
''Based on the steepness fragility index according to Eq. J66j. 
c Fit to Eq. H61|l for OTP was obtained by restricting log(ro) to 
be equal to —14. 
d From Ref. 
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Trom Ref. Ilia. The calculated dT g /dP refers to P = 300 MPa 
inst ead o f atmospheric pressure in case of OTP since the data in 
Ref. Ill9l app ly to high pressures only. 

-Trom Ref. Il20l . 

9 Using hig h-tem perature T, P-data for viscosity from Ref. Il2ll . 
h From Ref. H2I 

"Using the dielectric data at 0.1 MPa and 100 Mpa from Ref.ll23l. 



Combining Eqs. (J44J) , (|46|) . and (|59[) . the average relax- 
ation time of the hot region is then obtained as 



ln(r/r ) = 



2N max E(x) 



3T- {AE{x)/9)zC p (T) : 



(60) 



where -/V max is given by Eq. (|45[) and the molecular (mo- 
lar) conhgurational heat capacity Cp(T) appears in the 
denominator. For fragile liquids, E(x) is nearly indepen- 
dent of temperature down to the liquid-liquid transition 
point (x ~ 1, see Fig. [SJi) and can be considered as con- 
stant. This consideration yields the following equation 
for the relaxation time 



In (r/r ) = 



where the constant D is 



DT 



T - T'Cp(T) 1 



D = (9/4z)N c (N c - 1) 



(61) 



(62) 



and the temperature T" is described further below. 

The parameters D and T' in Eq. (|6ip can be considered 
as empirical fitting quantities for the sake of interpreting 
the experiment. According to Eq. (j3"12|) , the excess and 
configurational heat capacities of fragile liquids are equal 
to each other for fragile liquids. The 1/T scaling of the 
configurational heat capacity then results in the overall 
temperature dependence of the form 



ln(r/r ) 



DT'T 
T 2 - z\T> 



(63) 



This type of the temperature law was obtained for 
hard-sphere fluids by Jagla-^ who combined the Adam- 
Gibbs formalism with an empirical equation of state. 
Note, however, that the configurational entropy in Jagla's 
model has a 1 /T 2 temperature scaling inconsistent with 
the empirical 1/T law [Eqs. © and (|40|) ]. 



B. Application to experimental data 

Equations (jrTTj) and (|6"3"|) suggest some qualitative re- 
sults consistent with experimental observations. First, 
the model establishes a direct link between fragility 
and configurational/excess heat capacity. The configu- 
rational heat capacity in the denominator of Eq. (|6Tj) 
decreases at high temperatures, e.g. above the melting 
temperature T m . Therefore, the relaxation kinetics will 
change from super- Arrhenius at low temperatures to Ar- 



rhenius at high temperatures The extent of fragile 
behavior is controlled by A in the denominator of Eq. 
(|63p which, according to our thermodynamic analysis, 
is large for fragile liquids, resulting in curved Arrhenius 
plots. Note that parameter D in Eqs. ([61]) and (f63]) is 
essentially constant across all liquids studied (Table HIT)) 
in contrast to parameter D in the VFT equation [Eq. ([2])] 
which correlates with fragility. Finally, Eq. (|6ip predicts 
a return to the Arrhenius behavior— on passing below T g 
because of the drop of Cp which defines T g . This feature 
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DT7(T-T'C C P X ) 



FIG. 11: Dielectric relaxation time vs T (upper panel) and vs 
DT' '/{T-CpT') (lower panel) with D and T' listed in Table 
IIHI In the upper panel, the solid lines are fits to Eq. (|61[) and 
the dashed lines, indistinguishable from the solid lines on the 
scale of the plot, are the VFT fits. The intercept gives the 
physical value of 10~ 13 s. 
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FIG. 12: Adam-Gibbs plot of dielectric relaxation time vs 
(TS C {T))~ . The configurational entropy is calculated from 
the 1G model with the parameters listed in Table [TTJ The 
points are experimental dielectric data from Ref. I27I . the solid 
lines are obtained by using VFT fits of the experimental data 
from Refs. [8a and |89. 



has previously been unique to the Adam-Gibbs equation. 

In addition to providing a qualitatively correct picture, 
Eq. (|6ip performs surprisingly well in fitting the di- 
electric relaxation time o 27 ' 88 ! 89 of both intermediate and 
fragile liquids (Fig. fTTI and Table ITlT|) . as well as of simula- 
tion data (see below). Fits of experimental dielectric re- 
laxation times to Eq. (jBTj) , with to, D, and T" considered 
as fitting parameters, are indistinguishable from VFT fits 
in all cases studied (Fig. [TTJ upper panel). ln(r) plotted 
against DT'/(T - T'Cp x ) also yields straight lines and 
physically reasonable intercepts (Fig. 1111 lower panel). 
The quality of the linear correlations is as good as for 
the Adam-Gibbs plot (cf. Figs. fTTI and [T2|) and in some 
cases is even better (MTHF). In all fits, experimental 
Cp x (T) were used instead of Cp{T) suggested by the de- 
nominator of Eq. (j6"Tj) . In the 1G model, Cp x ~ Cp for 
fragile liquids. 

The nominator parameter D in Eq. ()61|) can be related 
to the average number of excitations in the hot region via 
Eq. UB2J. From TableHTTJ we obtain AT C ~ 10 - 20. With 
the usual molecular diameter of a ~ 5.5 A, this number 
projects into a length-scale of ~ 1.5 nm for a spheri- 
cal cluster or larger if relaxation is facilitated through 
chains of molecules^ This length is in general accord 
with current estimates of spatial heterogeneities in super- 
cooled liquids.*^ The activation barrier scales quadrati- 
cally with the number of units in the hot region compared 
to the linear scaling with the size of the cooperatively re- 
arranging region of the Adam-Gibbs theory. However, 
there is no divergent length-scale in the present formula- 
tion. 

Activated dynamics considered here do not anticipate 
the divergence of the relaxation time at any finite tem- 
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FIG. 13: Experimental Kauzmann temperature Tj^ p plotted 
against the temperature To alc at which the relaxation time in 
Eq. (|6ip diverges. This temperature is calculated by solving 
Eq. (|64[) in which temperature T' is taken from Tabic [TTT1 



perature, and the divergent solution given by Eqs. (|6ip 
and ([55)1 is an artifact of the approximate integration 
used in Eq. ([59]) . However, the mathematical solution of 
the equation 

To = T'C C P (T ) (64) 

can be associated with the VFT temperature To often 
reported experimentally. The temperature To calculated 
from this equation (Fig. [13)1 is almost equal to the ex- 
perimental Kauzmann temperature T^ p for all liquids 
listed in Table iLHl Equation (|61[) is therefore consistent 
with the empirically documented accord between To and 

rrnOXp 43 

1 K — 

For fragile liquids, Eq. (JH3J) simplifies to [Eq. (J32J), 
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CI 



C% 



r 



T 2 




TgC° p (Tg) 



(65) 



The temperature T" calculated from Eq. (|65p . with the 
thermodynamic parameters from Table HT1 is indeed close 
to the direct fit of Eq. ([61]) to experimental relaxation 
data (Table 1111)1 . The parameter T" can also be related 
to the steepness fragility index [Eq. I|33[)]. Assuming T _1 
scaling for Cp(T) in Eq. ([BTjl . one gets 



T' 



Cp(T g ) m + m„ 



(66) 



where 



log(r(T g )/r ) 



16. Estimates of V 



based on this equation are also listed in Table HTT1 Since, 
for fragile liquids, T" is related to observable quantities 
through Eq. (|65p . one can derive the following relation 
for the kinetic fragility 



m 1 + (T /T 5 ) 5 



1-(W) 2 ' 



(67) 



This equation, which was previously derived by Ruocco 
et alr^ holds reasonably well for fragile liquids (m > 
53). It makes a connection between the steepness index 
and the ratio Tq/T s recommended by some authors as a 
measure of fragility!^ 

Since the configurational heat capacity of fragile liquids 
is linearly related to the configurational entropy s c = 
Sq — cp(T), Eq. (|6ip can be re- written in terms of the 
configurational entropy 

ln(r/r ) 



DT' 



T + T'(S c -zs c ) 



(68) 



The current formulation still provides a link between the 
relaxation time and configurational entropy, but the alge- 
bra is different from the Adam-Gibbs relation [Eq. ((SJ)]. 
The validity of either functional dependence of the re- 
laxation time is often tested by calculating the slope of 
the glass transition temperature with pressure. Equa- 
tion (f6"Tj) suggests that T/T' — Cp(T) remains constant 
at the glass transition temperatures measured at differ- 
ent pressures. Following Goldstein ) 16 i 44 the condition 
d(T/T'-Cp{T)) = provides the slope of the glass tran- 
sition temperature with pressure, dT g /dP. One needs to 
calculate the derivative of T" over pressure and the pres- 
sure and temperature derivatives of Cp. The pressure 
derivative gives (dCp/dP)T = —T g V g Aap/kB, where V g 
is the molecular volume of the glass. The difference of 
the isobaric expansivities of the liquid and glass, Aap, is 
in the range ~ 5 x 

10 _ 4 K -i (Ref s . |i|[iJJ) 

allowing one 

to neglect the pressure derivative of Cp. One then gets 
^ = ^^(l-T'(dCp/dT) P )- 1 . (69) 



Equation (|69[) should be compared to what follows from 
the Adam-Gibbs scaling^ 
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FIG. 14: The liquid-liquid coexistence line ("LL-line") and 
the line of maximum heat capacity ("Widom line"^^) in the 
(T, so) parameters plane (a). Thermodynamic fragile behav- 
ior is limited by the critical point (T c ,sq c ) [Eq. (|37[) ], The 
vertical arrow shows the cooling path along which the relax- 
ation time is calculated using Eq. (|61|l (b). T max denotes the 
temperature at which the of heat capacity (dash-dotted line 
in (b)) passes through a maximum, C*p ax = CV(T max ). All 
the plots have been generated at constant en, A and z. 



Most 



liquids^ 5 * 122 



dT g /dP = T g V g Aa p /(S c + Cp). 



(70) 



data, available for strong/intermediate 

— i— indicate that Eq. (|7D|) adequately 

describ es the exp eriment. Since dTo/dP > 

(Refs. lmHl22lll27f ). one can expect from Eq. (giD 

that dT'/dP > 0. The actual numbers for this 

derivative can be obtained from T, P relaxation 
data^i^^^^^ The resmts Qf thege ca ,_ 

culation are given in Table IIII1 The relaxation times 
were fitted to Eq. (fSTj) at ambient pressure with D and 
T 1 considered as fitting parameters and experimental 
Cp x (T) and ambient pressure. The parameter D was 
then kept constant at elevated pressures and dT' jdP 
was evaluated from the fit. The results of these calcula- 
tions are in reasonable agreement with reported dT g /dP 
and reproduce the drop of this derivative in going from 
fragile to intermediate glass-formers ( Table [TTT]) . 



IV. DISCUSSION 

The model developed here describes the thermody- 
namic properties of glass-formers in terms of statistics of 
excitations from a single-energy state of the ideal glass 
to a Gaussian manifold of configurationally unjammcd 
states with higher energy and entropy. The model sug- 
gests that the thermodynamic signatures of fragile liq- 
uids, in particular a sharply increasing configurational 
heat capacity close to the glass transition, can be identi- 
fied with the existence of a thermodynamic phase tran- 
sition ("LL-line" in Fig. PHa). usually hidden below the 
glass transition temperature. 12 i 81 i 130 The strong/fragile 
behavior is distinguished within the model by two param- 
eters, the trapping energy parameter A and the entropy 
gain for a single configurational excitation sq, the ratio of 
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them equal to the experimental Kauzmann temperature 
[Eq. (|38|) ]. Both parameters significantly decrease when 
going from fragile to intermediate/strong liquids, with 
the parameter A showing the most significant change (Ta- 
ble |TT| . Fragile liquids arc therefore characterized by a 
broad range of configurationally excited states, whereas 
the energy distribution for intermediate/strong liquids is 
very narrow (Fig. Q]). Because of low values of either A 
or so (or both of them) strong/intermediate liquids fall 
in the range of temperatures and excitation entropies 
[T > T c , so < s 0c Eq. (j37|) ] not allowing a first-order 
transition. 
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FIG. 15: Illustration of the temperature behavior of the av- 
erage basin energy (4>) and the distribution of basin energies 
P(4>) in non-fragile and fragile liquids. The width of the basin 
energy distribution scales as 1/ \fN and requires small systems 
to be observed in simulations. 



The critical point (T c , so c ), separating fragile from non- 
fragile liquids fFig. [T^k). makes the descent into the en- 
ergy landscape qualitatively different for them. The av- 
erage basin energy of non-fragile liquids first starts to 
drop from the high-energy plateau according to the 1/T 
law and then inflects into an exponential decay. The 
distribution of basin energies is narrow, and its maxi- 
mum shifts to lower energies with cooling (Fig. [TSb.b). 
This behavior, often observed in simulations of binary LJ 
liquids ; 30 ' 33 ' 68 ' 95 is well reproduced by the present model 
which places these fluids into the strong/intermediate 
category As an example, we show in Fig. \W\ the fit 
of the 1G model to Sastry's data on binary LJ mixture 
(BL JM) i 28 i 131 The calculations were done for the highest 
density studied in simulations at which the BLJM liquid 
is most fragile. Equations (f2"B")) and (f3"Tj) were used, re- 
spectively, for the configurational entropy and heat ca- 
pacity, while Eq. (|61[) was applied to the simulated diffu- 
sion coefficient. The excitation energy (eo = 219 K) and 
the trapping energy (A = 71 K) obtained from the fit pro- 
duce the following critical parameters: T c = 35.5 K and 
sqc = 4.2 (accidentally, T c is close to T^ p ). On the other 
hand, the excitation entropy from the fit, Sq = 0.32, and 
the temperature range T > T c together put the BLJM 
liquid into the category of intermediate fragility. With 
such low excitation entropy, the population of the ex- 
cited state is close to 0.5 at high temperatures. From 
Eq. (22]) one then gets (z = 1) S Q = sg/2 + ln(2) ~ 0.85. 
The same number, 0.85, is reported by Sastry for the top 
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FIG. 16: Configurational entropy (a), configurational 
(constant- volume) heat capacity (b), and diffusion coefficient 
(c) of the 80(A) :20(B) BLJM liquid at pa\ A = 1.35, where 
<r aa is the LJ diameter of the component with higher con- 
centration. The thick solid lines are the results of simulations 
from Ref. [H arbitrari ly lim ited to T > 100 K. The diffusion 
coefficients from Ref. 1 13 ll (circles, c) are fitted to Eq. (|6ip 
(solid line) and to the VFT relation (dashed line). The two 
fits are almost indistinguishable on the scale of the plot. The 
fitting to 1G model [Eqs. (J2SJ, GO, and (JHTJ] results in the 
following fitting parameters: eo — 219 K, A = 71 K, sfj = 0.32, 
T' — 92.8 K, and D — 3.1. The thermodynamic Kauzmann 
temperature Tk is zero when calculated in the 1G model (the 
thin lines in a and b). 



of the enumeration function s^ ax i2£ 

The behavior of fragile liquids is quite different. The 
configurationally excited state is almost entirely popu- 
lated (x ~ 1) in the entire range of experimentally ac- 
cessible temperatures down to the glass transition. The 
average basin energy is almost temperature-independent, 
and the only effect of lowering the temperature on the 
distribution of basin energies is its narrowing accord- 
ing to the fluctuation-dissipation theorem (Fig. fT5b.d). 
The invariant population of the configurationally excited 
states is abruptly terminated at the equilibrium first- 
order transition or (on supercooling) at spinodal instabil- 
ity at which point it drops to a low value representative 
of an entropy-poor, low-temperature phase as experimen- 
tally confirmed for TPPf^ 

The fit of the model to experimental excess heat ca- 
pacities and entropies of molecular glass-formers (Ta- 
ble [TTJ) has pushed the liquid-liquid phase transition be- 
low the glass transition. For some substances,— ' 129 ' 134 ' 138 
the thermodynamic liquid-liquid transition can be found 
above the glass transition. As an example of this sit- 
uation, the 1G model is applied to excess entropy and 
heat capacity of laboratory water. 133 ' 135 The thick lines 
in Fig. [TTh.b show the results of measurements on bulk 
samples; 1 ^ 3 - while open circles refer to samples confined in 
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FIG. 17: (a) Excess entropy for bulk laboratory water 
("bulk"), nanoconfined water ("confined"), and configura- 
tional entropy for SPC/E water ("SPC/E"). The thick line 
and open points refer to bulk^ and nanoconfined exper- 
iments, respectively. The closed points refer to constant- 
volume simulations.- 1 ^ (b) Experimental excess heat capac- 
ity for bulk (thick line 1 ^) and hard nanoconfined 1 ^ (open 
points) water. The thin lines in (a) and (b) are fits to the 
thermodynamic 1G model [Eqs. (f27]l - (f29]) ] . The configura- 
tional heat capacity of SPC/E water in (b) is calculated from 
the 1G model based on the fit of the configurational entropy 
in (a), (c) Simulation (constant- volume) results for the Debye 
relaxation time (closed points) and diffusivity (open points) 
of SPC/E waterfi^ Diamonds show experimental diffusivity 
of laboratory water? 1 ^ The solid lines are fits to Eq. (|61[) 
with the configurational heat capacity obtained by numerical 
differentiation of S c for SPC/E water (see Table HVl for the fit- 
ting parameters) . The diffusivity curves have been arbitrarily 
shifted to fit to the scale of the plot. For all three panels, the 
closed points refer to simulations and open points are used to 
indicate the laboratory experiment. 



3 nm pores of silica gelJ^ The two sets of data coincide 
at high temperatures. Because of the limited tempera- 
ture range of the bulk water data the fit to the thermody- 
namic 1G model is done for nanoconfined water (Table 
IIV[) . The 1G model predicts a weak first-order phase 
transition at Tll = 229 K (P = 1 atm). This number 
is very close to previous estimates for the temperature 
of crossing the "Widom line" at which the heat capacity 
gets a maximumii^ and more recent measurements plac- 
ing this point at 223 K^. After the kink at T LL , the 
entropy slowly decays to zero at Tk — * 0. 

The behavior of excess thermodynamic parameters of 



TABLE IV: Fit of the 1G model to excess 
thermodynamic a 133 ' 135 and diffusivity 13 ' of laboratory 
water and to configurational entropy, diffusivity, and Debye 
relaxation time of SPC/E water from constant-volume 
simulation a 134 ' 136 ' 140 (p = 1.0 g/cm 3 ). The heat capacity 
and excess entropy of nanoconfined watep 1 ^ are used for the 
thermodynamic fitting. The energy parameters are in K, D 
and z are dimensionless. 

Liquid ep A sp s 0c z D T' log(r / s) 

Bulk/Confined water 1298 649 2.8 2.0 3.0 66 a 19 
SPC/E water 751 251 1.3 6 4.0 3.0 282 c 5.7 -13.4 

193 d 7.6 



"From fitting the diffusivity of laboratory water according to Price 
et alML 

''Configurational component of the excitation entropy. 
c From fitting the Debye relaxation time. 
d From fitting the diffusivity. 



the laboratory water is compared to configurational pa- 
rameters of SPC/E water obtained from constant-volume 
simulations of Starr et qZ, 134 . 140 The fit of Eqs. pg j) -(|29"l ) 
to the configurational entropy of SPC/E water (Fig. [T7a 
and Table [PV|) yields a low excitation entropy sg = 1.3, 
below the critical value of s 0c = 4.0 [Eq. ([37])]. The 
heat capacity curve then shows a broad pre-critical peak 
reminiscent of, but broader than, the heat capacity of 
nanoconfined water (Fig. [TTb). As in the case of the 
BLJM liquid, the excited state population, characteriz- 
ing the configurational manifold, is slightly above 0.5 at 
high temperatures resulting in the plateau of the config- 
urational entropy at (z = 3) So — .z(s§/2 + ln(2)) ~ 4. 
This number is close to the estimate for the top of the 
enumeration function given by Sciortino^ So ~ 3.7 at 
P ~ 0. 

The comparison of fits to laboratory and simulation 
data shows that the configurational part of the excita- 
tion entropy Sq makes about 46% of so for SPC/E water 
and only about 30% for o-terphenyl (Fig. [5]). These com- 
parable fractions of configurational components suggest 
that it is the change in the vibrational density of states 
(around the boson peak at ~ 30 cm -1 ) caused by con- 
figurational excitations that is primarily responsible for 
the fragile behavior of glass-formers^i 

The configurational entropy from simulations can next 
be used to fit the relaxation times by applying Eq. (|61[) . 
This procedure is straightforward for the Debye relax- 
ation time and diffusivity obtained from constant-volume 
simulations^^ at the same conditions as the data shown 
in Fig.fTTk (see Table HVl for the fitting parameters). The 
application to laboratory data is less obvious because of 
possible differences between configurational heat capac- 
ities at constant volume and constant pressure. Never- 
theless, diffusivity of laboratory water-i^l can be fitted 
to Eq. (|rjTj) , although with the fitting parameters (Table 
IIVP requiring more fragile behavior (cf. triangles to open 
circles in Fig. [TTb). 

Equation (fSTj) is the central result of our modeling of 
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relaxation. It is also consistent with the empirically docu- 
mented accord between the VFT and Kauzmann temper- 
atures fFig. U3p . The model combines the thermodynam- 
ics of configurational excitations with the notion that as- 
sembling excited units in clusters is required for activated 
relaxation* 1 - 4 ^ The thermodynamics of glass-formers is 
thus determined by the statistics of excitations, while 
dynamics probe more rare events of clustering of excita- 
tions. One of the signatures of thermodynamic driving 
force behind relaxation is the realization of a dynamic 
fragilc-to-strong transition close to the point where the 
heat capacity in the denominator of Eq. (|6ip reaches its 
maximum, Cp ax . This sub-critical maximum is indeed 
predicted by the thermodynamic 1G model is the vicin- 
ity of the critical point (T c ,sq c ) (Fig. [14a). When the 
cooling path goes close to the critical point, the frag- 
ile non-Arrhenius dynamics changes to strong Arrhcnius 
dynamics after passing the line of the heat capacity max- 
imum ( "Widom line" r^i Fig. [Mb) . A similar behavior 
has recently been observed for nanoconfined wate n 139 i 142 
Experimental evidence supports a link between the vis- 
cosity and configurational thermodynamic o 28 ' 33 ' 34 ' 143 but 
this does not necessarily mean that activated events are 
driven by the entropy. Equation (|61[) still connects the 
relaxation time to configurational entropy [see Eq. (|68|) ]. 
but is based on standard arguments of activated kinet- 
ics in terms of excess kinetic energy accumulated at a 
"mobile unit" . The model thus puts focus on the kinetic 
energy and its fluctuations as the main driving force of 
liquid relaxation.— This notion brings relaxation of su- 
percooled liquid in general accord with approaches de- 
veloped to describe activated events in chemistry where 
accumulation of sufficient kinetic energy along a reaction 
coordinate, and not the entropy, is viewed as the gen- 



eral mechanism behind activated transitions. Since re- 
laxation is understood in terms of clusters of excitations, 
the model can potentially bridge to more macroscopic ar- 
guments of elastic theories considering a region of shear 
displacements shoving the environment from a relaxing 
unit4£ 

The present model does not anticipate a divergent 
length-scale and also avoids the need to consider the mo- 
saic interface energies of Ref.H^, which recent simulations 
have failed to support It remains to be seen whether 
the formalism of configurational excitationS ) 53 i 55 ' 56 i 57 ' 144 
"dressed" with the relevant excitation thermodynamics, 
provides a sufficient basis for the theoretical modeling 
of supercooled liquids. While the model yields a set of 
simple equations for some basic properties commonly re- 
ported for glass- formers [e.g., Eqs. (|32|) . (|38|) . (|37|) . and 
(|oTj) ] , it should be recognized that the mean-field charac- 
ter of the model may push it too much into a strong first- 
order type of thermodynamic transition. Some "soften- 
ing" of the model including fluctuations around the mean 
field would result in greater flexibility, but this would 
only come at the expense of an increase in the number of 
parameters. 
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